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ABSTRACT 


The goal of this thesis is to investigate the behavior 
of a thermoacoustic prime mover filled with argon. Data was 
collected for three different operating conditions: 1) prior 
to the onset of acoustic oscillation, 2) immediately after 
onset and 3) after the prime mover had reached steady state, 
at mean gas pressures of approximately 80 kPa, 205 kPa and 305 
kPa. Measured values of 1/Q are compared to predictions of a 
theory based on Swift's work (Swift, G. W., "Thermoacoustic 
Phones," Jj. Acoust. Soc. Am. Vol. 84, October 1988). After 
onset, Q is determined from measurement of the rise time of 
the oscillations. The data set collected is the first to span 
onset and results show good agreement with theory. It is 
concluded that the theory can be used to provide reasonable 


input for a model to predict the steady state waveforms. 
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[us INTRODUCTION 

A thermoacoustic prime mover is a heat engine; DE 
converts stored energy into useful work, in the form of sound. 
The basic operation and construction of a prime mover can be 
explained with reference to Fig. 1. For our purposes, a prime 
mover can be considered to be a closed resonator comprised of 
five sections. A temperature gradient is applied across the 
prime mover stack via the ambient and hot heat exchangers. 
(The terms "ambient" and "hot" refer to the temperatures of 
the heat exchangers and the portions of the resonator adjacent 
to them.) The (prime mover) stack and heat exchangers are 
parallel plates spaced by a few thermal penetration depths and 
Oriented such that the plates lie along the axis of the 
resonator. The hot and ambient sections are circular tubes -- 
Open at the heat exchanger end and closed at the other end. 
The net absorption coefficient of the prime mover is the sum 
of the thermal and viscous wall losses in each of the five 
Sections. Attenuation in the gas is negligible. 

As the temperature difference imposed across the stack 
increases from zero, the thermal losses in the stack decrease 
and eventually become negative, representing gain. At a 
sufficiently high temperature difference, the gain of the 
stack overcomes the remaining losses and the prime mover 


spontaneously goes into self oscillation. This phenomenon is 
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Figure 1. Major components of a prime mover. (Not to 
scale) 


known as onset. Once onset is reached, the acoustic amplitude 
of the resonator grows rapidly, typically reaching a few 
percent of the mean pressure. This discussion leads to the 
conclusion that three distinct regions of prime mover behavior 
are of interest: 1) below onset, 2) above onset during the 
buildup of oscillations to steady state, and 3) steady state. 
The purpose of this thesis is to investigate each of these 
regions in an argon-filled prime mover. 

This thesis 1S one in a sequence that have investigated 
thermoacoustic prime movers. The goal of this research is to 
understand the behavior of prime movers in each of the three 
regions of operation. In 1989, Lin began a quantitative study 
of a helium-filled prime mover. He investigated the region 
below onset as well as documenting the nature of steady state 
oscillations. Lin found that the steady state oscillations in 
our prime mover are highly nonlinear, exhibiting a rich 
harmonic structure. (Lin, 1989) Therefore, it became clear 
that, to completely understand the steady state region, it is 
also necessary to understand the generation of harmonics and 
energy dissipation in high amplitude standing waves. In 1991, 
Chen began such an investigation in a simple standing wave 
tube (Chen, 1991). Also in 1991, Bowers investigated the 
initial buildup of oscillations above onset (Bowers, 1991). 

As stated earlier, the purpose of this thesis is to 
investigate the three regions of behavior in an argon-filled 


prime mover. Ones reasons tor Enis. as ~ to, adds eto our 


understanding of prime movers by using a gas other than 
helium. Also, Lin and Bowers used slightly different prime 
movers so direct comparison of their data is not possible. We 
Will use the same prime mover configuration below and above 
onset so that we can have a complete data set from zero 
temperature difference through onset and beyond. 

Following a discussion of applicable theory, the 
experimental setup and procedures used to collect data will be 
presented. The results and their analysis follow, along with 
a discussion of sources of error. Finally, a summary of this 
experiment and its conclusions will be made along with 


recommendations for subsequent research in this area. 


II. THEORY 


Theoretically, the performance of a prime mover can be 
investigated through the calculation of its acoustic power 
output. However, below onset and during the initial buildup 
of oscillations, it is easier to measure the quality factor of 
the prime mover. 

ne met aceustie power ourput ZW. can besrelated to the 
quality factor, Q. Below onset, Q is easily determined from 
the prime mover’s frequency response. Just above onset, Q can 
be determined from the rise time of the oscillations. This 
Section contains a discussion of the theory necessary to 
calculate 1/Q in both regions. This discussion is a summary 
of that found in Bowers (Bowers, 1991, pp. 10-15) which draws 
heavily from Swift (Swift, 1988). 

The reciprocal of the quality factor, 1/0, of the prime 


mover can be expressed as 


EAS un | (1) 
Q OL or 





pere w is the angular frequency and E, the stored energy of 
the resonant tube. W is determined by integrating the power 
per unit volume, w(x), over the volume of the prime mover, as 


follows, 
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S(x) is the cross-sectional area of the prime mover at x. 


In the stack, w(x) is 
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P, and k are the peak pressure amplitude and wave number =a. 


the standing wave, whereas a and 0 represent the thermodynamic 


sound speed and Prandtl number of the working fluid. J =a] 
Also, VT, is the mean temperature gradient along the wee a 


plates and T is the mean temperature of the plate anda: 








function of x. f, is defined as 
ር ከ) 
camp] e 
f = E (4) 
1 + 7) - 
( 7) 5 
for parallel plates. ár the viscous penetration depuis 
defined as 
See eee (5) 
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u is the fluid’s dynamic viscosity. A similar relationship 
exists between f, 6 and y,, although 6 is defined as 


AE 2K | (6) 
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Reais the thermal conductivity of the working fluid, 
ep cdctincs the relative neat Capacity or tne gas compared 


to the plates 
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where p, and Cc, are the mean density and isobaric specific 
Meat or the working fluid and p., c. and é represent similar 
qualities for the stack plates. 1 is one-half the stack plate 
thickness. In the heat exchangers, w(x) is identical to 
Huston (3), except AT — Q0. In the ambient and hot ends, 
HEUS Given by Equation (3) with ATO and | = 0 (there 
are no plates). Also, the functions f, and f, take on 
different forms involving the Bessel functions. 

E. ls determined by integrating the time averaged energy 


ST 


density over each prime mover section, 
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where S is the cross-sectional area of the resonator and q is 
the phase of the standing wave as it enters each section. 

These equations use the properties of the working ክስክስ ን 
and the initial conditions of the prime mover to predict 1/9 
as a function of VT. acpossatucgs toes 

The resonance frequency is calculated by neglecting the 
dependence of the sound speed, a, on dissipation. In other 
Words a is assumed to be the thermodynamic sound speed. Next, 
we require that acoustic pressure antinodes exist at each end 


of the prime mover which amounts to requiring that 


fi (x) dx - nr. (9) 


k depends on x through k(T (x)). 


III. EXPERIMENT 


A. EXPERIMENTAL SETUP 
The prime mover and temperature and pressure 
monitoring/control systems are same ones used in Bowers' 
thesis. Only pertinent details will be given here. 
1. The Prime Mover 
As discussed earlier and shown in Fig. 1, the prime 
mover is composed of five sections. 
a. The Hot End 
The hot end of the prime mover consists of a nickel 
tube, 0.63 cm thick, 5.5 cm long and having an inner diameter 
of 3.82 cm. A nickel cap is welded into the top of this tube 
to provide a rigid boundary. To monitor the hot end 
temperature, a small hole is drilled into the end cap and 
fitted with an access tube for a K-type thermocouple. The 
lower end of this assembly holds the hot heat exchanger. 
b. The Hot Heat Exchanger 
The hot heat exchanger is also composed of nickel 
and like the hot end, has an inner diameter of 3.82 cm. 25 
nickel plates, 0.045 cm thick and 0.762 cm long are mounted 
along the heat exchanger’s longitudinal axis. The spacing 
between these plates is 0.104 cm. 304 stainless steel spacers 


mounted to one side of each plate help maintain this spacing 
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and prevent deformation under high heat conditions. The 
hot end/heat exchanger assembly 1S welded to the prime mover 
stack. 
c. The Prime Mover Stack 

The prime mover stack is housed in a thin walled 
304 stainless steel tube, 3.5 cm long, with an inner diameter 
of 3.82 cm and a wall thickness of 0.05 cm. The stack 
consists of 35, 304 stainless steel plates, 0.028 cm thick, 
Spaced 0.077 cm apart. As in the hot heat exchanger, the 
Spacing between these plates iS maintained by using 304 
stainless steel spacers mounted to one side of each plate. 

d. The Ambient Heat Exchanger 

A constant ambient temperature is maintained at the 
bottom of the stack by using 2 identical copper heat 
exchangers. Each 1S 1.02 cm long and consists of 25 copper 
plates, 0.025 cm thick and spaced 0.105 cm apart. The inner 
diameter of each heat exchanger is 3.82 cm. The two heat 
exchangers are spaced 0.15 cm apart. The top heat exchanger 
is mounted in contact with the lower end of the prime mover 
and positioned so that its plates form SM with the 
stack plates. This maximizes fluid flow out of the stack by 
ensuring that the heat exchanger plates do not block the gaps 
between the prime mover plates. The second heat exchanger is 
not necessary for the purposes of this thesis. It will be 


used in later investigations. 
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A constant working temperature is maintained within 
the heat exchangers by mounting them ina brass water jacket. 
This water jacket also contains a small hole used for 
pressurizing the prime mover. The hot end heat 
exchanger/stack assembly is soldered to the top of the water 
jacket. 

e. The Ambient End 

The ambient end of the prime mover consists of a 
3.82 cm ID, 88.11 cm long copper tube, flanged at the Bose 
to accommodate an end cap. Two different end caps are used. 
Below onset the end cap houses an electret driver and an 
electret microphone, so that the frequency response of the 
tube can be determined. The end cap used to collect data 
after onset holds only a microphone to record the amplitude of 
the standing wave. 

A butterfly valve is located 38.0 cm from the 
bottom of the ambient heat exchanger. Closing the butterfly 
valve effectively shortens the length of the prime mover and 
requires a greater differential temperature to drive the prime 
mover to onset. With the butterfly valve closed, a 
differential temperature across the stack sufficient for onset 
can be established without actually generating a standing wave 


(Bowers, 1991, p. 18). 
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2. Control and Monitoring Systems 
a. Temperature 

The hot end of the prime mover is fitted with an 
Omega Engineering HBA Model 202040 heater. Power is supplied 
to the heater by a General Radio Type 100-Q variable 
transformer which controls the heater voltage. The heater 
coils are wrapped tightly around the hot end of the prime 
mover and stainless steel shim stock is used to maintain good 
thermal contact. A Neslab RTE-110 constant temperature bath 
supplies cooling water to the water jacket and to plastic 
tubing which encircles the entire length of the ambient end of 
the prime mover. 

Temperatures along the length of the prime mover are 
monitored by a Keithley Model 740 System Scanning Thermometer. 
A K-type thermocouple probe is inserted into the access hole 
in the hot end until that it makes contact with the hot heat 
exchanger. Three E-type thermocouple junctions are affixed to 
the outside of the ambient end of the prime mover with 
cyanoacrylate adhesive. These junctions are placed below the 
ambient heat exchanger, below the butterfly valve and above 
the pressure transducer flange at the bottom of the ambient 


end. 
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b. Pressure 

The prime mover in this thesis uses argon gas as 
the working fluid: The gas pressure is controlled by a 
Matheson 8H-580 pressure regulator with a working range of 
O - 1400 kPa. Prior to pressurization, the prime mover is 
evacuated using a Welch Model 1400 Duo-Seal Vacuum Pump for 
several hours to ensure all residual gases are removed. 
Pressure within the prime mover is monitored via the pressure 
sensing line in the waterjacket assembly. An external tubing 
and valve network allows the pressure sensing line to be 
isolated from the gas supply cylinder and vacuum pump. The 
pressure is monitored by an Omega PX304-300AV pressure 
transducer having a sensitivity of 3PSIA/MVDC. The output of 
the pressure transducer is read by a Hewlett Packard 3457A 


digital multimeter. 


B. EXPERIMENTAL PROCEDURE 
As power is supplied to the hot end of the prime mover, 
the differential temperature, AT , increases from zero Eca 
level sufficient to bring about onset. Data was acquired in 
the three regions of prime mover operation discussed in the 
Totroduction: 
1. Below Onset 
a. Data Acquisition 
A schematic diagram of the apparatus used to 


measure frequency response of the prime mover below onset is 
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Shown in Fig. 2. Data acquisition is controlled using GPIB 
hardware and a program written in Microsoft QuickBasic and 
executed on a Standard 286 PC-AT. After entering the expected 
resonance frequency, desired bandwidth, number of frequency 
increments and starting system pressure, the PC directs a 
Hewlett-Packard 3324A Synthesized Frequency Generator to 
output a 10 Vpp sine wave beginning at the lower end of the 
desired bandwidth. This signal is sent to the reference input 
of a Stanford Research Systems 530 Lock-In Amplifier, a 
Kikusul COS6100A Oscilloscope and a Techron 7520 Power Supply 
Amplifier which amplifies the signal and sends it to the 
driver. 

The driver is a locally constructed electret type 
which sets up a standing wave in the prime mover. An electret 
microphone is used to sense the amplitude of the standing 
wave. The output of the microphone is first directed through 
a locally constructed preamplifier which boosts the signal by 
a factor of 100. The amplified signal is then sent to the 
lock-in amplifier and oscilloscope. 

The oscilloscope provides a visual display of the 
reference signal and the received signal, while the lock-in 
amplifier filters the received signal in a narrow bandwidth 
around the reference frequency. It sends the in-phase and 
quadrature components of the signal to the HP 3457A 
multimeter. These values are read by the PC, via a GPIB 


interface, and the amplitude of the signal is determined. 
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Figure 2. Diagram of Apparatus Used to Collect Data Prior 
to Onset 
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After the acoustic amplitude data and the 
temperature of each thermocouple are written to a worksheet, 
as shown in Fig. 3, the frequency is incremented. This 
procedure is repeated until the end of the bandwidth is 
reached, at which time the ending system pressure is recorded. 
Data was recorded for both the fundamental and second modes of 
the prime mover at the three mean gas pressures of interest. 

b. Data Processing 

Q is obtained by applying the same fitting function 

Lin used to determine the resonant response of a simple tube. 


nar function is 


GOD) | (9) 
kL 





me 
Dre ape 
de) 


where P is the normalized amplitude of the acoustic pressure 





in the tube and A is a scaling constant. The subscript "o" 
indicates a measured resonant value. KL is defined as 
kL = i| Y + 7 3 (40) 
Wo (20) 





By performing a least squares fit of P versus f, values for 


f A and Q can be determined. (ham. 1989 po. 228 Data is 


o! 
processed using macro routines in Microsoft's Excel program. 
Figure 4 shows a typical plot resulting from the fitting 


procedure. 


py 


ለ | C 
File Name: 31030011.XLS 
3 |Microphone Gain: 100 

F. Gen. Amplitude: 10000mV 
5 [Center Frequency: 161Hz__ | | | | . 
_6 |Start Pressure: 306.7159 KPa_ | 0 | .: | . .: | . 
7[:.:...  መዕጩሬልመዉ፡ሜጫጮጫሜማ፡ማ፡፡ስ፡ሜይ፡፡፡ወ፡፡ ጋ. 
515 [A | ኅበ O 
. 099483) 1472206116... 1... oo ስሕ 

rn ቤር: 
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E 
|12| 160.0006| 10.72061| 0.391125| 3089... 206] 204) 204 
|13| 160.2008, 13.37106| 0.365257 | 3089, 206 20.4) 204 


14] 160.4005 15.48943| 0.434295 308.9 206 | 204 

|15| 160.6005. 19.19547| 0.565027 308.9) — 206] 20.4) 20 

|16| 160.8008 23.49876 0.827011 308.9 20.6 20.4 20.4 

|17| 161.6009| 27.43512| 1.046905 4 3089, 206] - 204. - 204 

| 18| 161.2008 29.3778, 1.366104 3089 206 _ 204 204 
161.4003. 27.86031| 1.690933 308.9 206 204 

|20| 161.6009 25.73375| 1.989192 3089 206 

21| 161.8006! 22.04756 2.247514 308.7 20.6 20.4 20.4 

22| 162.0009) 17.27475| 2.421878 3086. 206| 204) - 20.41 

1231: o spas $|we — | ee | tc c | 

| 24 |End Pressure: 297.1531 KPa .‹ጠው | 


25 | Average Pressure: 301.9345 KPa, | | | 

26 B cou EM a ስክ 
Pavg — |Ttop — |Sdev  |Tbott  |Sdev |Fo Jo | 
28 | 301.9345] 308.8364| 0.102691] ^ 20.6| 426E-07| 161.26 147.76 





Figure 3. sample worksheet used to record data prior to 
onset. 
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Figure 4. Sample plot of data obtained prior to onset. 


To compensate for any loss of pressure during the 
run, the mean gas pressure is obtained by averaging the 
starting and ending system pressures. The mean hot and 
ambient end temperatures are the averages of the temperatures 
recorded at each incremental frequency. 

2. Determining Onset 
Unlike helium, which exhibits a well-defined onset 


temperature when used as the working fluid, it was more 
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difficult to identify the experimental onset temperature in 
argon. 

As onset was approached, a very low amplitude acoustic 
wave was detected by the sensing equipment. The amplitude of 
this wave increased slowly, requiring, in some cases, several 
minutes to build to steady state. The exact reason for this 
behavior 1S unknown. The temperature differences required for 
onset are shown in Table 1. These values were determined by 
slowly increasing AT until onset occurred, then observe Ale 
behavior of the prime mover at several temperatures around 


this reference point. AT is the mean value of the 


Onset 


temperatures for which onset was observed. At onset, 1/Q= 0. 





ATE) Bio EC cM OS ES 


Table 17) SUMMARY OF ONSET TEMPERATURES 


3. Above Onset 
a. Data Acquisition 
A schematic diagram of the apparatus used to 
collect data above onset is shown in Fig. 5. The pressure 
transducer used for this phase of data collection is an 
ENDEVCO Model 8510B-5 piezoresistive microphone with a linear 
response range of + 5 PSI. Data is collected above onset 


using the procedures outlined in Bowers’ thesis. With the 
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E oure 5: Diagram of apparatus used to collect data after 
onset. 
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is established and allowed to stabilize. When the butterfly 
valve is opened, acoustic oscillations begin to build up in 
the prime mover. The resulting acoustic wave is sensed by the 
microphone and routed through a locally constructed, 
differential amplifier having a gain of 10. This signal is 
sent to a Nicolet 310 digital oscilloscope for recording and 
storage. (Bowers, 1991, p 20) The oscilloscope is triggered 
to begin recording data immediately after onset occurs, 
ensuring that the acoustic signal is small. This minimizes 
the presence of harmonics in the recorded signal and ensures 
that the temperature distribution in the stack 1s unaltered by 
Significant acoustic heat transport. 
b. Data Processing 
The Q of the prime mover above onset is determined 


by the relationship, 
Ww 
C" 13 
Q 2B Jos, 


where o is the angular resonance frequency and ß is the 
temporal absorption coefficient. To ascertain the value of 
these variables, the data collected above onset is digitally 
processed using the techniques developed by Bowers (Bowers, 
1991 5 DD 222-277) 

The binary signal stored by the Nicolet 310 is 
converted to ASCII code to be processed using Mathworks MATLAB 


program. Initially, a discrete Fourier transform, DFT, is 
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performed and the positions of the resonance frequency and any 
harmonics are determined. A raw waveform is shown in Fig. 6, 
while Fig. 7 shows the DFT of this data. 

The original signal is multiplied by a complex 
exponential term which effectively shifts the fundamental 
Dreguenc vto DC (i.e. 0 Hz). zA 5th order, 30 Hz Butterworth 
filter is then applied to the shifted data set to isolate the 
fundamental frequency. Figure 8 shows a plot of the shifted, 
filtered data revealing that only the envelope of the original 
data remains. As shown in Bowers' thesis, the end effects of 
the Butterworth filter can be avoided by disregarding the 
first 300 and last 300 points of filtered data. The 
exponential envelope of the original waveform is linearized by 
taking the logarithm of the normalized amplitude, so that 
aleast mean squares fit can be performed. The slope 
determined by this fit corresponds to the desired quantity, f. 

4. Steady State 
a. Data Acquisition 

Collection of data once the prime mover has reached 
steady state uses the same experimental SES as was used to 
collect data above onset. With the butterfly valve open and 
17 Werjejenkt AT, established Losınaucee self oscillation che 

prime mover is allowed to reach equilibrium and the resulting 


waveform is recorded by the Nicolet oscilloscope. 
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Figure 8. Plot of shifted, filtered data collected during 
onset. 





In order to prevent aliasing, the digital sampling 


rate must conform to the Nyquist criteria, 


fe > 2 (12) 


Nyquist max * 


where f is the minimum sampling frequency and f,. is the 


Nyquist 
highest frequency of interest in the signal being sampled. 
(Strumm, 1988, p 54) Since the fundamental frequency of the 
prime mover 1S approximately 160 Hz and there will be at most 
Menar monicsidetectaplelin cheisteady state = comanda 
E uyauist ^ 3520 H2. This corresponds to a minimum sampling rate 


of 308 us per point. We used 200 us per point. 
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Since we are interested in the absolute pressure amplitude 
of the harmonics of the steady state waveform, a known test 
Signal is also recorded, using the same oscilloscope settings. 
The DFT of this signal will provide a reference when analyzing 
the original data. 

b. Data Processing 
Using Mathworks MATLAB program, a DFT of the 
reference signal and the recorded waveform are performed, as 
shown in Figs 9 and 10. In the interest of clarity, Fig. 9 
shows only the first 1000 points of the original signal's 4000 
point DFT, while only the first half of the DFT of the 


recorded waveform is shown in Fig. 10. 
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Figure 9. Plot of DFT of reference signal used for 
analysis of steady state data. 
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Figure 10. Sample plot of DFT of data collected during 
Steady state. 
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IV. RESULTS 


Data were taken at 80 kPa, 205 kPa and 305 kPa for A MEE 
ranging from O'C to approximately 425'C, the self-imposed safe 
Working limit for the prime mover assembly. Below onset, data 
for the fundamental and second modes were collected at all 
three pressures. However, due to temperature constraints, 
only the two lower pressures were investigated above onset and 
at steady state. 

1. Results of 1/Q Measurements 

Figures 11 and 12 show the combined results of the 
data collected below onset and immediately above onset for 
pressures of 80 kPa and 205 kPa. Fig. 13 shows the rece 
obtained below onset at 305 kPa. For discussion purposes, it 
should be pointed out that 1/Q is proportional to the net 
attenuation in the prime mover. 

The first point of discussion is that the data in 
Figs. 11 and 12 show a smooth transition from below onset, up 
to and beyond onset. This feature indicates that the methods 
of determining 1/Q in each region are consistent. These data 
are the first to show the progression from below onset to 
above onset. It should be pointed out that the above onset 


data were obtained at low acoustic amplitudes. This ensures 
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Figure 11. Plot of combined 1/Q vs Al data obtained- berow 
and immediately after onset. P = 80kPa. 
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ን ከርበ] 1 Plot of combined 1/0 vs AT data obtained below 
oand immediately after onset. P = 205 kPa. 


29 


200 
Delta T ( C) 





Figure 13. Plot of 1/0 vs AT, for data obtained befowNcm c 
in the fundamental mode. P = 305 kPa. 


that the effects of acoustic heat transport can be ignored, 
just as they were below onset. The results for the second 
mode are shown in Figs. 14, 15 and 16. All of these data were 
obtained below onset. 

The agreement between data and theory is good, 
however, there is a tendency to underpredict the attenuation. 
This trend is especially evident in the second mode. One 
possible explanation for this is the fact that the model does 
not include the presence of the butterfly valve. The 
magnitude of this effect should be roughly proportional to the 
additional surface area introduced by the valve. The valve is 


a circular disk of negligible thickness with a diameter of 
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Pigure 14. Plov of 1/Q vs AT for data obtained below onset 
in the second mode. P - 80 kPa. 
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Ingure 15. Plot of 170. vs AT. for data obtained below onset 
in the second mode. P = 205 kPa. 
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Figure 16. Plot of 1/0 vs VT, for data obtained below onset 
in the second mode. P = 305 kPa. 


approximately 38 cm. The total SUBÉSCS area is, therefore, 
approximately 24 cm. The total surface area of the 
primemover, including heat exchangers and prime mover stack, 
is approximately 2300 cm. Therefore, the valve e 
approximately 1 percent more surface area to the prime mover. 
Although a more detailed calculation of the effect of the 
valve is needed, it is unlikely that it alone can account for 
all the error. It should be pointed out that the same trend 
is evident in helium data. (Atchley, 1991) 

Another possible explanation concerns the fact that 
and T 


i are recorded at the top of the hot heat exchanger 


hot cold 


3,2 


and the bottom of the ambient heat exchanger, respectively. 
These temperatures may not accurately reflect the actual 
temperature difference across the stack. This suggests that 
the actual Ar ulee than that measured, an error that would 
in effect shift the data to the left, providing a much better 
overall agreement. The magnitude of this error is difficult 
to measure. Measurements performed to try and estimate the 
error revealed that the results were very sensitive to the 
degree of contact between the heat exchangers and the stack. 
A discrepancy of + 10°C is certainly not unreasonable. 
However, without knowing the degree of contact in the prime 
mover, any estimates of the magnitude of the error are 
speculation. 

Another point of comparison has to do with the 
dependence of the difference between theory and experiment on 
An Figures 11 and 12, and, to a lesser extent, Fig. 13, 
e ate the agreement 156 worse at higher values of AT. 
This trend is absent in Figs. 14, 15 and 16. A possible 
explanation for this trend may the irregularities in stack 
geometry. The theory assumes plate spacing is uniform, 
however, this is not the case. One would expect this effect 
to be most pronounced when ፅ is comparable to y,. Recalling 
eae definition of & from Egn. (6); tnis would be true for low 
pressures and frequencies. This follows exactly the trend 


noted. 
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2. Other Sources of Error 
In addition to the specific sources of error 
previously mentioned, there are others associated with the 
apparatus and procedures which apply to all data collected. 
a. Monitoring Systems 


The measurements of mean gas pressure, P and 


iit 
system temperatures are subject to instrumentation errors and 
errors due to averaging. T. andi m are usedi to determine AT 
while T. and P, are inputs to the equation used for theorem i 
prediction of prime mover response. 

For example, when investigating the response of the 
fundamental mode at P ካ9ሻ፡ ቭ EUM AN SON the average 
value of Ta OVer the course of data collection is 172.2 + 
0.5°9@ and T. “= 20:2 4 02103) C mas arest Tc ee 
0.5 C; an error of 0:55 percent. p, is the average cio 
mean system pressures at the start and finish of each data 
collection run. Actual system pressure varies over the course 
of collecting data due to small leaks in the prime mover. 
Typically, the system lost one percent of pressure at 205 kPa, 
two percent at 305 kPa and gained 0.2 percent at 75 kPa. 

The prediction of prime mover behavior for a given 
pressure and) mode requires መሪ np These 
parameters are obtained by averaging their values for an 


entire data set. For example, in the fundamental mode at 
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kPa, differences of 1.5 and 1.0 percent, respectively. 


Furthermore, the sensors used to monitor these 
parameters are subject to design errors. According to the 
manufacturer’s specifications, the error 11) K-type 


thermocouples is + 5°C while the E-type thermocouples are 
subject to an error of + 2°C. The Omega pressure transducer 
has an error of + 0.5 percent, based on manufacturer's data. 
b. Determining the Resonance Frequency 
Prior to onset, the fitting function used to 
Jeremie ss the resonance frequency,  f. is accurate Tto 
approximately 0.1 percent. Above onset, f. ais determined by 
analyzing the DFT of the sampled waveform. እ sampling rate of 
5 usec for a 4000 point record corresponds to a frequency 
range of 2000 Hz. This results in an accuracy of 0.5 Hz per 
PONE OF 073 percent tor {= 160°H2. A Jarger error in the 
pxcuUeney occurs as a result of “calculating f- It was stated 
in Chapter II that the computation neglects dispersion due to 
attenuation. The result is that the theory slightly 
ሠ ተ However, the: maximum error jt approximately 
2 percent, 
3. Steady State 
Figures 17 through 20 show the waveform and DFT of 
data collected after the prime mover had reached steady state 


above onset at approximately 80 kPa. Figs Il A 129517) 
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correspond to a AT, of 382.8 C, while Figs 19 ando 
correspond to a AT, of 419.8'C. An extensive set of data is 
presented in the Appendix. It is seen from these figures that 
the waveforms become highly nonlinear at higher temperature 
differences. These results are consistent with Lin'/s data 
obtained in helium. A future goal of this research is to 


predict these waveforms. 
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Figure 17. Plot of waveform obtained at steady state. 
P = 80.1 KPa. AT, = 382.8 “Er 
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Figure 18. Plot of DFT of data obtained at steady state. 
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Figure 19. Plot of waveform obtained at steady state. 
P = 80.3 kPa. AT, - 419.8 °C. 
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Figure 20. Plot of DFT of data obtained at steady state. 
P = 80.3 kPa. VT, = 419.8 'C. 


38 


V. CONCLUSIONS 


The purpose of this thesis was to investigate the behavior 
of an argon-filled thermoacoustic prime mover. The reasons 
for this investigation were two-fold. First, the data 
collected below and immediately above thermoacoustic onset 
allowed a verification of theoretical predictions begun in 
previous research. Since the previous work used helium as the 
working fluid, this experiment varied the input parameters by 
using argon, instead. Additionally, this data is the first to 
Span onset. Our results indicate that the theory provides a 
good agreement with measurement, although some discrepancies 
do exist. At lower mean gas pressures, the prime mover’s 
behavior begins to deviate from theory at higher values of 
ቹፒ.) 111071 may belque to imperífections in the prime 
mover stack geometry. There is also a tendency to 
underpredict the attenuation. 

Data was also obtained at steady state. These data 
show that the waveforms are highly nonlinear. One of the 
ultimate goals of our research is to be able to predict the 
shape of these waveforms. At this time we are working with 
a model that predicts the waveforms of high amplitude standing 
waves in simple tubes. The inputs to the model are the 
resonance frequency and the Q of the modes of the tube. Once 


a prime mover has gone into oscillation, it is no longer 
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possible to measure these parameters. We will have to rely on 
theoretical values. 

The good agreement between the measured and predicted 
values of 1/Q indicate that we will be able to provide a 
reasonably good input for the model. This discussion leads to 
a somewhat subtle, but important, point. Below onset, Q is 
determined by fitting the measured frequency response to the 
ideal frequency response of a simple tube. The quality of the 
agreement indicates that, to a good approximation, the prime 
mover behaves like a simple tube having a uniform cross 


section and temperature distribution! 
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APPENDIX 


PLOTS OF STEADY STATE DATA 
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Figure 21. Plot of waveform obtained at steady state. 
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Figure 22. Plot of DFT of data obtained at steady state. 
P = 81.1 kPa. AT, = 398.3 'C. 


42 


q 
E 
2 
E 
= 
< 
© 
። 
= 
E 
« 


0.015 0.02 





Time (sec) 
Figure 23. Plot of waveform obtained at steady state. 
P = 81.2 kPa. AT, = 394.2 'C. 
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Figure 24. Plot of DFT of data obtained at steady state. 
Be UD 2 KPa: AT 0394.9 "€. 
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Figure 25. Plot of waveform obtained at steady state. 
P = 80.3 kPa. AT, = 411.4 °C. 
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Figure 26. Plot of waveform obtained at steady state. 
P = 80.3 kPa. AT = 411.4 “Cc, 
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Figure 27. Plot of waveform obtained at steady state. 
P = 206.7 kPa. AT, = 384.7 °C. 
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Figure 28. Plot of DFT of data obtained at steady state. 
P = 206.7 kPa. AT, = 384.7 'C. 
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Figure 29. Plot of waveform obtained at steady state. 
P = 206.5 kPa. AT = 395.1 Os 
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Figure 30. Plot of DFT of data obtained at steady state. 
P = 206.5 kPa. AT, = 395.1 °C. 
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Figure 31. Plot of waveform obtained at steady state. 
P = 206.7 kPa. AT, = 398.7 'C. 
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Figure 32. Plot of DFT of data obtained at steady state. 
P= 206.7 kPa. AT. — 398.7 'C. 
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Figure 33. Plot of waveform obtained at steady state. 
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Figure 34. Plot of DFT of data obtained at steady state. 
P = 204.8 kPa. NT = ንታ Gr 
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Figure 35. Plot of waveform obtained at steady state. 
DEL Ou Dada = dc Dc. 


Amplitude (Arbitrary Units) 


ሪኣ 


100 200 300 400 500 600 700 800 900 





Frequency (Hz) 


Figure 36. Plot of DFT of data obtained at steady state. 
pu 207.) Da Aq 306.0. 9 
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